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ABSTRACT 

We present submillimetre and mid-infrared imaging observations of five fields centred on 
quasi-stellar objects (QSOs) at 1.7 < z < 2.8. All 5 QSOs were detected previously at sub- 
millimetre wavelengths. At 850 (450) ftm we detect 17 (11) submillimetre galaxies (SMGs) 
in addition to the QSOs. The total area mapped at 850 iim is 28 arcmin^ down to RMS 
noise levels of 1 — 2 mJy/beam, depending on the field. Integral number counts are computed 
from the 850 ptm data using the same analytical techniques adopted by 'blank-field' submil- 
limetre surveys. We find that the 'QSO-field' counts show a clear excess over the blank-field 
counts at deboosted flux densities of 2 — 4 mJy; at higher flux densities the counts are 
consistent with the blank-field counts. Robust mid-infrared counterparts are identified for all 
four submillimetre detected QSOs and ~ 60 per cent of the SMGs. The mid-infrared colours 
of the QSOs ai-e similar- to those of the local ULIRG/AGN Mrk 231 if placed at 1 < z < 3 
whilst most of the SMGs have colours very similar to those of the local ULIRG Arp 220 at 
1 < z < 3. Mid-infrared diagnostics therefore find no strong evidence that the SMGs host 
buried AGN although we cannot rule out such a possibility. Taken together our results sug- 
gest that the QSOs sit in regions of the early universe which are undergoing an enhanced 
level of major star-formation activity, and should evolve to become similarly dense regions 
containing massive galaxies at the present epoch. Finally, we find evidence that the level of 
star-formation activity in individual galaxies appears to be lower around the QSOs than it is 
around more powerful radio-loud AGN at higher redshifts. 

Key words: galaxies: formation - galaxies: starburst - galaxies: evolution - galaxies: high 
redshift - submillimetre - QSO: general 



1 INTRODUCTION 

The popular hierarchical model of galaxy formation predicts that 
elliptical galaxies found today in the cores of rich clusters formed 
at high redshifts and at rare high-density peaks of the dark matter 
distribution. Within these regions, gas-rich proto-galaxies merge 
together rapidly, forming stars at a high rate. The same reservoir 
of gas used to build the stellar mass can also fuel the growth of 
the supermassive black holes (SMBH) found dormant in the centre 
of the galaxies at low redshift (e.g. Kauffmann & Haehnelt 2000; 
Page et al. 2001; Hopkins et al. 2008). Since massive star forma- 
tion is known to be a dusty phenomenon, the optical light from 
many young galaxies should be highly obscured. Dust enshrouded 
star formation is, however, a luminous phenomenon in the millime- 
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tre through far-infrared wavebands where the star-light absorbed 
by dust grains is re-emitted. Indeed, in the last decade, deep sur- 
veys at submillimetre (e.g. Small, Ivison & Blain 1997; Hughes et 
al. 1998; Bales et al. 1999; Scott et al. 2002; Coppin et al. 2006) 
and millimetre (e.g. Greve et al. 2004; Bertoldi et al. 2007) wave- 
lengths have identified a high-redshift population of massive, dusty 
galaxies that are undergoing extreme starbursts (e.g. Smail et al. 
2002). Likewise, hotter dust heated in the circumnuclear environ- 
ment of an AGN will emit at mid-infrared wavelengths. Therefore, 
submillimetre and mid-infrared observations can be combined to 
study activity in forming galaxies due to dust obscured starbursts 
and AGN. 

To-date there has been no published survey at (sub)millimetre 
or far-infrared wavelengths of a randomly selected region of sky 
over an area sufficiently large to ensure that very rare structures in 
the early universe are contained within it. Although large surveys of 
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this type will be available soon from instruments such as SCUBA-2 
on the JCMT and SPIRE/PACS on-board the Herschel Space Ob- 
servatory, published surveys have so far been limited to an area of 
~ 0.5 deg"^ (Coppin et al. 2006; Austermann et al. 2010; WeiB 
et al. 2009). One of the largest published submillimetre surveys is 
the SCUBA HAlf-Degree Extragalactic Survey (SHADES; Mortier 
et al. 2005; Coppin et al. 2006) which, due to technical problems, 
mapped only 0.2 deg^ discovering ^ 100 submillimetre galaxies 
(SMGs). Whilst this survey provided one of the largest and most 
uniformly selected samples of SMGs ever assembled, the area of 
sky surveyed is extremely small by most standards. For example, 
ACDM simulations predict that a dark matter halo of mass greater 
than a few times 10^* at z; = 2 will evolve to become a very 
rich cluster of galaxies by 2 = (e.g. Farrah et al. 2006). However, 
since the surface density of such massive haloes at z > 2 is only 
of order 0.01 deg^^ (Evrard et al. 2002), a wide-area survey would 
be required in order to contain one by chance. 

One promising method of locating rare structures is to tar- 
get the fields of high-redshift active galactic nuclei (AGN). Given 
their huge luminosities, such objects must contain a SMBH even 
at very high redshifts, and therefore must represent some of the 
most massive objects in existence at their epochs. They should thus 
act as signposts to rare high-density regions of the early universe, 
and can be observed efficiently with existing technology. Several 
groups have exploited this technique to search for proto-clusters 
forming around high-redshift radio galaxies (e.g. Kurk et al. 2000; 
Pentericci et al. 2000; Ivison et al. 2000; Stevens et al. 2003; Small 
et al. 2003; De Breuck et al. 2004; Greve et al. 2007; Venemans et 
al. 2007) and quasi-stellar objects (QSOs; Stevens et al. 2004; Prid- 
dey, Ivison & Isaak 2008). Most of these studies have concentrated 
on extreme objects, those with the highest luminosities and at the 
highest redshifts, and all have found over-densities of star-forming 
galaxies (typically Lya emitters or SMGs) in the fields of the AGN. 
Here we extend this work to investigate the environments of typical 
QSOs at submillimetre and mid-infrared wavelengths. The QSOs 
are selected over the redshift range 1.7 < z < 2.8 where the 
star-formation rate density and accretion luminosity density of the 
universe peak. This work is motivated by the desire to study the 
coupled formation of stellar mass and black holes in galaxies that 
should evolve to become typical ellipticals in moderately rich clus- 
ters by the present day. Since the bulk of the field SMG population 
are also found at these redshifts (Chapman et al. 2003, 2005) we 
can also compare the properties of galaxies found in different envi- 
ronments. 

This paper presents submillimetre and mid-infrared imaging 
data for 5 QSO fields. We extract submillimetre catalogues and in- 
tegral number counts, and compare these with results from blank 
field surveys. We then characterise the mid-infrared properties of 
the QSOs and the SMGs discovered in their fields. 

A Hubble constant Hq = 70 kms~^ Mpc~^ and density pa- 
rameters SIa ~ 0.7 and Q,n-i = 0.3 are assumed throughout. 



2 SAMPLE SELECTION 

Our targets are drawn from the sample of 28 hard X-ray selected 
QSOs described by Page, Mittaz & Carrera (2001). The objects 
have X-ray spectra harder than the X-ray background, consistent 
with moderate absorption (20.5 < log(NH/cm^^) < 23) and dis- 
play broad optical emission lines. In previous work (Page et al. 
2001; Stevens et al. 2005) we have observed 20 of these objects us- 
ing SCUBA in photometry mode. The remaining 8 objects were not 



targetted with SCUBA due to observing constraints so this sample 
should be representative of the population. Our sample for subse- 
quent mapping comprises those QSOs with photometry mode de- 
tections of > 5 mjy/beam. This flux density threshold is based 
on limitations imposed by current sensitivity limits and by confu- 
sion noise, but also means that any source detected at 2: > 1 must 
be either an ultraluminous or hyperluminous far-infrared galaxy if 
the emission is from dust, and if the temperature of this dust is 
typical of such objects detected at low redshift. The observations 
of RX J094144.51 + 385434.8 were published by Stevens et al. 
(2004) but are included here for completeness. The 5 QSOs are 
listed in Table [T] along with some basic parameters. They are se- 
lected to (1) span the redshift range 1.7 < z < 2.8 and (2) have 
0.5 — 2 keV X-ray luminosities within ~ 0.7 dex of Lx, the break 
in the X-ray luminosity function. 



3 OBSERVATIONS AND DATA REDUCTION 
3.1 SCUBA observations 

SCUBA jiggle-map observations at 450 and 850 /im were made 
at the JCMT between 2003 July and 2004 December. These ob- 
servations are summarized in Table [2] We used a 64-point jiggle 
pattern to fully sample the image plane at both wavelengths but em- 
ployed a number of different chop/nod strategies. For RX J163308, 
RX J094144 and RX J121803 we chopped and nodded East- West 
with throws of 30, 30 and 45 arcsecs respectively. For these objects 
the off-positions fall at a fixed position on the sky; a real source 
thus has negative images of itself to the East and West, separated by 
twice the chop throw but with half the amplitude. For RX J005734 
and RX J124913 we chopped and nodded in azimuth by 30 arcsec 
so the off-positions rotate with the field throughout the observation. 

Each mapping observation, taking typically one hour to com- 
plete, was bracketed by pointing observations on bright point-like 
objects. Pointing offsets were in most cases < 2 arcsec and in all 
cases < 3 arcsec. Focus was checked throughout the night but par- 
ticularly around sunset and sunrise when rapid changes are known 
to occur. 

Atmospheric transmission was determined with skydips and 
checked against the JCMT water vapour radiometer. Maps of 
Uranus, Mars and/or CRL618 were used for flux density calibra- 
tion while maps of bright blazars were used to measure the beam 
profile. At 850-/im the calibration was found to be stable enough 
to apply a single calibration to all observations; our combined flux 
calibration factor (FCF) was 221 ± 15 Jy/V/beam giving a cali- 
bration uncertainty of about 7 per cent. At 450-/im the FCF varies 
with the atmospheric conditions because the filter profile contains a 
relatively strong water line, but also depends quite strongly on the 
thermal state of the antenna (see Jenness et al. 2002). At 450 ^m 
we thus applied our best estimate of the FCF to each observation 
in turn; FCFs varied from 220 to 490 Jy/V/beam. The calibration 
uncertainty at 450 ^m is ~ 20 per cent. 

Data were reduced with standard STARLINK SURF routines up 
to the rebinning stage (extinction correction, flat-fielding and de- 
spiking). For each mapping observation the individual bolometer 
timestreams were inspected and noisy integrations/noisy bolome- 
ters were blanked. We then used adapted versions of the SURF tasks 
SETBOLWT and REBIN to make bolometer noise weighted mean 
signal and noise maps on a 2 arcsec pixel scale (see Stevens et al. 
2004 for more details). These maps were then convolved with a 
Gaussian of full width half maximum 4 pixels (850 fim) or 3 pixels 
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Table 1. The sample of QSOs mapped with SCUBA. Coordinates are those of the optical QSO. Absorption corrected X-ray luminosities are given for the 
0.5 — 2 keV band whilst far-infrared luminosities are calculated from submillimetre photometric observations at 850 /xm (see Stevens et al. 2005 for more 
details on how these luminosities were calculated). 
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Table 2. Observing log of the SCUBA mapping observations. For each field 
we give the date of observation, the number of 64-pt jiggle maps where 
10 maps is equivalent to ~ 20 mins of integration, and the opacity range 
(Twvm) as given at 225 GHz from measurements with the JCMT water 
vapour radiometer. Observations of RXJ094144. 51 + 385434.8 were pub- 
lished by Stevens et al. (2004). 
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(450 /im) giving final resolutions of 14.2 (850 /xm) and 8.6 arcsecs 
(450 /itn). The convolved maps with signal-to-noise (S/N) contours 
are shown in Fig.[T] 



dithers with frame times of 100 or 200s depending on the predicted 
brightness of the target. Frame times for the MIPS observations 
were between 3 and 10 s. These observations are summarised in 
Tabled 

The IRAC frames were processed by the Spitzer Science Cen- 
ter using the standard pipeline version S 14.0.0. The images were 
then cropped to only include the region of constant exposure time; 
this allows catalogues to be constructed from regions with simi- 
lar noise properties. Source extraction was performed with SEx- 
TRACTOR (Berlin & Arnouts 1996) with a fixed aperture of 3.5" 
and a local background fit. Aperture corrections of 1.6 and 2.13 
were applied to the 4.5 and 8.0 /im flux densities respectively, in 
accord with the values determined from the Spitzer First Look Sur- 
vey (Lacy et al. 2005). 

For MIPS observations, the 24 /im BCD images were fitted 
with a plane after first masking out the bright objects. This plane 
was then subtracted from each BCD frame resulting in a flat back- 
ground image. These images were then combined with the IRAF 
task IMCOMBINE using the astrometric solution from each indi- 
vidual BCD image. This procedure resulted in a much improved 
24 fim image, free of the jail-bar pattern common to the PBCD 
images supplied from the Spitzer pipeline. Again, only the central 
3x3 arcmin^ region of the images was used to detect sources. 
We used SEXTRACTOR to find sources and the PHOT script in IDL 
to measure their flux density within a 13 arcsec radius aperture. All 
flux densities were then subjected to an aperture correction of 1.167 
(see Seymour et al. 2007). 



4 THE SUBMILLIMETRE CATALOGUES 

Our philosophy is to match our analysis as closely as possible to 
that performed by the SHADES consortium (Coppin et al. 2006) 
thus allowing a direct comparison of the number counts determined 
for blank fields and those in the highly biased regions of the Uni- 
verse surrounding our QSOs. Most of the analysis in this section 
pertains to the 850 /im data. For the 450 fim data we extract a cat- 
alogue at a fixed S/N threshold but do not compute number counts. 



3.2 Spitzer observations 

We observed each field with the IRAC (at 4.5 and 8.0 ^m) and 
MIPS (at 24 fim) cameras on-board the Spitzer Space Telescope. 
IRAC observations were made between 2005 March and 2005 July 
while the MIPS observations were made between 2005 March and 
2005 December. For the IRAC observations we used 5-pt Gaussian 



4.1 Source extraction 

We adopt the source extraction algorithm described by Scott et 
al. (2002) which performs a maximum-likelihood fit to the flux 
densities of all potentially significant peaks identified on the con- 
volved maps. For each source we took the convolved S/N map and 
searched for all peaks with S/N> 1.5. A peak is found if a 10 x 10 
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Table 3. Observing log of the Spitzer observations. For each field we give the date of observation and the on-source integration time (Tint) f™ each instrument. 
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arcsec grid centred on any pixel fails to find another pixel in the 
grid with higher S/N. In practice, this means that real sources sep- 
arated by less than a beam width can be missed by the extraction 
algorithm. This is the case for RX J094144 where we fail to re- 
cover two sources which are clearly real because we see them in the 
higher resolution 450 /.im data but are blended together at 850 /im. 
We then took the unconvolved signal map and convolved it with 
appropriately scaled beam maps centred at these locations. The am- 
plitudes of these beam maps were then adjusted until the overall 
of the fit was minimised. For two sources with azimuth chop/nod 
schemes the beam was reconstructed by going back to the observ- 
ing logs and calculating the position angle of the chop throw on 
the sky at the mid-point of each observation. We then added these 
together, weighted by the integration time for each observation. 

This method is suitable for extracting sources in fields where 
confusion is likely to be a problem; it also has the advantage of us- 
ing the off-positions to distinguish between real sources and noise 
spikes, and increases the S/N by considering the flux density in 
the off-positions. The method was used for all four independent 
SHADES reductions although note that 'Reduction A' of Coppin et 
al. (2006) searched for peaks on the convolved signal maps rather 
than the convolved S/N maps. We tried both methods but found that 
peaks towards the map edges were better recovered from the S/N 
maps. 

Uncertainties in submillimetre positions measured on SCUBA 
maps have been discussed extensively by Ivison et al. (2007) who 
show that the 1-a error circle is given by 0.9 6'a(SNR2 - S)-'''^ 
where 9a is the FWHM of the beam and SNR is the raw map signal- 
to-noise ratio. A typical source with SNR~ 3 — 4 at 850 /im would 
thus have a l-cr positional accuracy of 3 — 5 arcsec. 

4.2 Flux density deboosting 

A well known phenomenon common to all deep submillimetre sur- 
veys is that of flux density boosting. Real sources falling in regions 
of positive noise will have their intrinsic flux densities boosted, 
raising the likelihood of inclusion in the catalogue. This effect be- 
comes increasingly significant as the intrinsic flux density of the 
source approaches the limit of the survey. If not corrected for, flux 
density boosting will therefore populate the survey catalogue with 
an excess of low flux-density sources. 

We use the Bayesian flux density deboosting recipe of Coppin 
et al. (2005, 2006) which calculates a posterior flux density prob- 
ability distribution for each source assuming a prior distribution of 
real flux densities constructed from existing data. The prior is con- 
structed using artificial skies which are sampled with the same ob- 
serving technique used for the real data. Since we used two differ- 
ent chop/nod schemes we employ two priors; a two position chop 
scheme for the 3 fields observed with East- West chop throws and 
the SHADES observing pattern for the two fields that were chopped 



in azimuth. SHADES employed a 6 position chop/nod scheme, 
chopping in both right ascension and declination with chop throws 
of 30, 44 and 68 arcsecs. We note that while this method is not a 
true match to an azimuth chop scheme, where the chop position 
angle rotates on the sky throughout the observation, it is neverthe- 
less a better approximation than a two chop scheme. In any case, 
running the deboosting algorithm with either simulated observing 
technique resulted in the same catalogue give or take I or 2 sources 
close to the selection threshold. Deboosted flux densities are listed 
in TableEl 



4.3 Catalogue membership 

As a first cut we took all sources found by the extraction algorithm 
with a S/N> 2.0. This catalogue was then trimmed to include only 
those sources which have less than 5 per cent of their posterior 
probability distribution below mjy. This cut is exactly that used 
by Coppin et al. (2006) for the SHADES catalogue. In Fig.|2]we 
plot extracted 850-/im S/N against extracted 850-/im flux density. 
Those sources falling above the solid line, and that remain there af- 
ter deboosting, meet the above criterion and are included in the final 
catalogue (TableQ. In the 5 fields we find 21 sources at 850 fim 
including 4 of the 5 QSOs. At 450 /im we simply list those sources 
extracted with S/N> 3 (Tabled and do not attempt to deboost their 
flux densities. We assume that an 850 /im source is also detected 
at 450 ^m if their extracted positions are within 10 arcsec. Using 
this criterion we recover 8 of the sources detected at 850 in- 
cluding 3 of the QSOs. A further 6 sources detected at 450 /im are 
not detected at 850 ^m. Two of these objects lie in the RX J094144 
field and were missed at 850 /im due to blending as discussed in 
SectionO We note that the QSO RX J121803.82 + 470854.6 is 
not formally detected at either 850 or 450 ^m even though it was 
detected by Page et al. (2001) using SCUBA in photometry mode. 
This QSO is discussed further in Appendix A. 



4.4 Catalogue completeness 

Extensive simulations were performed to determine the complete- 
ness as a function of flux density. For each field we added false 
sources to the signal maps by scaling the normalised beam map to 
have an amplitude of 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 14, 16 and 18 mJy. 
We then attempted to recover each source with the extraction algo- 
rithm and repeated this procedure 1000 times per field. We consider 
a source detected if its recovered flux density is within a factor of 
2 of its input flux density, if its recovered position is within 7 arc- 
secs of its input position and if it makes it into our catalogue as 
described in the previous section. These criteria are those adopted 
by the SHADES survey. We plot the fractional completeness for 
each field as a function of input flux density in Fig. |3] The fits to 
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Figure 1. SCUBAmaps of (a) RX J005734.87-272827.4, (b) RXJ094144.5H-3854.34.8, (c) RXJ121803.82+470854.6, (d)RXJ124913.86-055906.2 
and (e) RX J163308.57 + 570256.7 at 850 (left) and 450 //m (right). The grey-scales show flux density while the contours, which show signal-to-noise, are 
calculated using the noise maps (see text). Contour levels are 2, 3, 4, 5 and 6 a except for RX J124913.86 — 055906.2 where they are 2, 3, 4, 5, 6, 8, 10 and 
12 a. These contoitfs refer to the significance levels on the pre-extracted maps as shown. Note that real sources will have higher significance when extracted 
with the beam profile because emission in the off-positons is also factored into the total signal-to-noise. 




Figure 1. - continued 



the data have the form Completeness — (S'")/(6 + cS") where 
S is the 850-/im flux density and a, b, c are free parameters. It is 
apparent that the completeness is a function of both the average 
noise level of the maps and the structure that they contain. We note 
that the fractional completeness converges to ~ 0.85 — 0.9 rather 
than 1 at high flux densities. This is due to the failure of the extrac- 
tion algorithm to recover input sources which fall very close to the 
field edges. The effect will be worse towards the top and bottom 
centres of the fields where both off-positions fall outside the area 
used for extraction. We confirmed that this is the case by perform- 
ing an 18 mJy simulation over the centre region of one of the maps. 
This simulation returned a fractional completeness of 1.0. How- 
ever, since real sources may also fall on the map edges in exactly 
the same way that randomly generated false sources do, we do not 
correct for this edge effect. One other feature of these simulations 
worthy of remark is that the shape of the completeness curve ap- 
pears to be sensitive to the adopted chopping scheme. Those fields 
observed with an azimuthal chop/nod scheme have relatively high 
completeness at low flux densities, presumably due to an absence 
of the deep negative off-positions that are produced by chopping 
at a fixed position on the sky (most fields have a bright QSO at 
their centre), but they are relatively less complete at high flux den- 



sities where the absence of these off-positions becomes more of a 
hindrance. 



5 INTEGRAL NUMBER COUNTS 

Number counts at 850 /im were calculated from the deboosted flux 
densities after first making a correction for the completeness of the 
catalogue for each field using the curves shown in Fig. |3] The to- 
tal area of the survey is 28.4 arcmin^. We follow Coppin et al. 
(2006) so that our counts are calculated in the same manner as 
the SHADES counts. The method uses the deboosted flux density 
probability distributions as input for a modified boot- strapping sim- 
ulation which estimates the number counts and their errors. This 
method allows the flux density of an individual source to be shared 
appropriately between adjacent bins. We quote counts in bins of 
width 2 mJy (QSOs included; Table[6) and 1 mJy (QSOs removed; 
Table IT). These counts represent the mean value in each bin with 
the quoted error being the frequentist 68% confidence interval in 
each flux bin based on boot-strapped Monte Carlo simulations, in- 
cluding Poisson counting errors. While including the QSOs gives a 
biased representation of the over-density in these fields, this method 
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Table 4. The 850-//m catalogue. Coordinates give the extracted position of the source. Also listed are extracted flux density and S/N together with deeboosted 
flux density (see text for details). Quoted uncertainties on the raw flux densities include contributions from the calibration uncertainty and the extracted 
signal-to-noise values while those on the deboosted flux densities are purely statistical. 



Field 




ID# 


RA 
(J2000.0) 




Dec 
(J2000.0) 




'S'sso 
(mjy/beam) 


5850 (deb.) 
(mJy/beam) 


Map S/N 


Note 


RXJ005734.78 


- 272827.4 


1 


00'' 57™ 35^ 


.21 


—27° 


28' 


27 


.9 


8.2 ±1.2 


'•'^-o.g 


7.9 


QSO 






2 


00'' 57™ 37= 


.79 


—27° 


28' 


57 


.9 


5.9 ± 1.6 


4 3+1.7 


3.8 




RXJ094144.51 


+ 385434.8 


1 


09'' 41™ 44= 


.84 


+38° 


54 


nail 

3d 


.9 


12.7 ± 1.4 


12 4+1-' 


11.4 


QSO 






2 


09'' 41™ 46= 


.53 


+38° 


CO/ 

53 


^ All 

54 


.7 


7.4 ± 1.9 


c 7+2.0 


4.1 








3 


Ogh 4im 42= 


.80 


+38° 


55 


n All 

34 


.8 


5.5 ±1.7 


o j+2.1 


3.3 




RXJ121803.82 


+ 470854.6 


1 


12'' 18™ 05= 


.51 


+47° 


08' 


55 


.1 


4.3 ± 1.3 


o j + 1.3 


3.5 








2 


12'' 18™ 07= 


.66 


+47° 


08' 


45 


.1 


4.5 ±1.4 


2.8i;;« 


3.3 








3 


12'' 18™ 05= 


.53 


+47° 


08' 


nail 

3d 


.9 


3.6 ±1.2 


2 1+1-3 
^•^-1.3 


3.1 








4 


12'' 18™ 02= 


.76 


+47° 


08' 


n All 

24 


.9 


3.5 ± 1.2 




3.1 








5 


12'' 18™ 04= 


.74 


+47° 


08' 


59" 


.0 


3.6 ± 1.2 




3.0 




RXJ124913.86 


- 055906.2 


1 


12'' 49™ 13= 


.99 


-05° 


59' 


19" 


.3 


11.2 ± 1.1 


ii.oi":^ 


13.9 


QSO 






2 


12'' 49™ 15= 


.45 


-05° 


58' 


37" 


.3 


4.6 ± 1.1 


3.8il;l 


4.3 








3 


12'' 49'° 17= 


.06 


-05° 


59' 


33" 


.3 


3.1 ±0.9 


9 O+0.9 


3.5 








4 


12'' 49'° 12= 


.37 


-05° 


58' 


33" 


.5 


3.2 ±1.0 


n Q+10 
^••^-1.0 


3.4 








5 


12'' 49'° 11= 


.56 


-05° 


59' 


59" 


.3 


3.2 ± 1.0 


2 1+1-1 


3.2 




RXJ163308.57 


+ 570258.7 


1 


16'' 32™ 59= 


.99 


+57° 


03' 


18" 


.7 


10.7 ± 1.4 


10.4^1-0 


9.5 








2 


16'' 33™ 08= 


.82 


+57° 


02' 


58" 


.7 


7.1 ±0.9 


6 q+"-^ 


9.3 


QSO 






3 


16'' 33™ 12= 


.74 


+57° 


03' 


26" 


.7 


4.5 ± 1.0 


3.8li;« 


4.5 








4 


16'' 33"° 10= 


.78 


+57° 


02' 


10" 


.7 


4.0 ±1.1 




3.7 








5 


16'' 33"° 03= 


.18 


+57° 


02' 


58" 


.7 


4.0 ± 1.2 




3.5 








6 


16'' 33"° 12= 


.61 


+57° 


02' 


32" 


.0 


3.4 ±1.1 


2.3il;i 


3.4 





Table 5. The 450-/im catalogue. Coordinates give the extracted position of the source. Also listed are extracted flux density and S/N. Quoted uncertainties on 
the raw flux densities include contributions from the calibration uncertainty and the extracted signal-to-noise values. The last column gives the likely ID if the 
450 /um position falls within 10 arcsecs of either the QSO or one of the objects extracted at 850 ;um. 



Field 


ID# 


RA 
(J2000.0) 




Dec 
(J2000.0) 




5450 

(mJy/beam) 


Map S/N 


Note 


RXJ005734.78 - 272827.4 


1 


00'' 57™ 35= 


.08 


-27° 28' 


34" 


.1 


32.4 ± 10.2 


4.1 


QSO 




2 


00'' 57™ 33= 


.13 


-27° 28' 


40" 


.1 


25.9 ±9.6 


3.2 




RXJ094144.51 + 385434.8 


1 


09'' 41™ 45= 


.06 


+38° 54' 


36" 


.7 


44.7 ± 10.2 


9.1 


QSO 




2 


09'' 41™ 46= 


.73 


+38° 53' 


56" 


.7 


41.2 ± 10.0 


7.3 


850#2 




3 


ogh 4im 4QB 


.73 


+38° 54' 


14" 


.9 


33.5 ± 8.8 


5.9 






4 


ogh 4im 44s 


.16 


+38° 54' 


22" 


.7 


28.8 ± 7.8 


5.5 




RX J121803.82 + 470854.6 


1 


12'' 18'° 05= 


.51 


+47° 08' 


47" 


.1 


35.1 ± 8.9 


6.4 


850#1 




2 


12'' 18™ 00= 


.80 


+47° 08' 


34" 


.9 


31.9 ± 10.7 


3.7 






3 


12'' 18™ 05= 


.11 


+47° 08' 


35" 


.1 


30.4 ± 11.0 


3.3 


850#3 




4 


12'' 18™ 10= 


.02 


+47° 08' 


41" 


.1 


24.3 ±8.8 


3.3 






5 


12h 18™ 07= 


.08 


+47° 08' 


52" 


.9 


21.9 ±8.3 


3.1 


850#2 


RXJ124913.86 - 055906.2 


1 


12h 49™ 13= 


.97 


-05° 59' 


21" 


.3 


34.9 ±8.1 


8.7 


QSO 




2 


12h 49m lis 


.03 


-05° 58' 


37" 


.5 


21.6 ±7.7 


3.4 






3 


12h 49™ 11= 


.03 


-05° 59' 


57" 


.5 


33.7 ± 12.5 


3.2 


850#5 
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• RX J005734 

* EX J094144 
» EX J181803 
■ RX J124913 
ft EX J 163308 



5 10 15 

Observed 850— /im flux density (mJy) 

Figure 2. Extracted 850 /im flux density against extracted 850 fim signal- 
to-noise (uncorrected for flux density boosting). All sources extracted with 
S/N> 2 are plotted; those that fall above the solid line have > 95 per cent 
posterior probability that their flux density is > mJy. These sources make 
it into our catalogue. Dotted lines show noise levels of 1 and 2 mjy/beam. 



gives the total number of submillimetre detected galaxies. This ap- 
proach is also analogous to that adopted by previous workers on 
QSO over-densities (e.g. Priddey et al. 2008), although we stress 
that the counts are only appropriate for X-ray absorbed QSOs that 
are pre-selected to be submillimetre luminous. 

Quoted number counts exclude the contribution from sources 
3, 4 and 5 in the RX J121803 field. All three have deboosted flux 
densities close to 2 mJy where the completeness estimate from our 
simulations is very low (< 1 per cent) and therefore subject to 
a large uncertainty. The implication is that number counts in the 
lowest flux density bins will be underestimated. If we assume that 
the true completeness of this field at 2 mJy is 5 per cent (a rea- 
sonable value from inspection of Fig. O and that the three sources 
contribute entirely to the lowest flux density bin then the number 
counts in this bin would increase by 70 per cent. However, no cor- 
rection is applied to the quoted counts. 

The counts (QSOs removed) are plotted in Fig. |4] which also 
shows number counts determined for the SHADES survey (Cop- 
pin et al. 2006). We note that counts determined at 870 /xm from 
a slightly larger survey conducted with the Large APEX Bolome- 
ter Camera (LABOCA) are consistent with those measured for 
SHADES at flux densities below 3 mJy and are about a factor 
of 2 below the SHADES counts at higher flux densities (WeiC et 
al. 2009). The dashed line is not a fit to the data but is simply 
the best-fit Schechter function to the SHADES counts multiplied 
by 4. From this plot we can see that the observed over-density of 
sources around the QSOs arises almost entirely from galaxies with 
deboosted flux densities in the range 2 — 4 mJy. At higher flux 
densities the counts in the QSO fields are matched very well to 
those determined from SHADES data for the field. 



Table 6. Combined integral number counts for the 5 QSO fields including 
the QSOs. Bins are quoted at the lower flux density bound. The bins are of 
width 2 mJy. SHADES counts are taken from Coppin et al. (2006). 



Flux density 850 fim integral counts, N(> S) 
(mJy) (deg-2) 

This work SHADES 



2.0 




2506lfo^ 


4.0 


2074+1060 


844+11^ 


6.0 


953+507 


3621^1 


8.0 






10.0 


q7c; + 190 
■^'^-237 




12.0 







Table 7. Combined integral number counts for the 5 QSO fields excluding 
the QSOs. Bins are quoted at the lower flux density bound. The bins are of 
width 1 mJy. 



Flux density 850 fim integral counts, N{> S) 
(mJy) (deg-2) 
This work 



2.0 


10420+3'' 


3.0 


■^°'^-1211 


4.0 


1423lg4 


5.0 


6361308 


6.0 


040+205 


7.0 


997+156 
-167 


8.0 


-175+146 


9.0 


130+130 



6 MID-INFRARED COUNTERPARTS 

Likely mid-infrared counterparts to the submillimetre galaxies are 
given in Table [8] For each field we searched within fixed radii of 
the extracted submillimetre positions. These radii of 5 arcsec at 4.5 
and 8.0 fim, and 10 arcsec at 24 /im were chosen to approximately 
match the theoretical error circles on the submillimetre positions 
added in quadrature with the uncertainty on the mid-infrared posi- 
tions, which are significant at 24 ^m (see Ivison et al. 2007 for a full 
discussion). The corrected P statistic (Downes et al. 1986) is given 
for each source found within the search radii. This quantity is given 
by P = 1 - exp[-P.{l + ln(Pc/P.)}] where Pc = ■kt'^Nt, 
P* — vrr^TV, r-g is the search radius, A^t is the surface density of 
objects at the limiting flux density of the catalogue, A'^ is the sur- 
face density of objects brighter than 5* where S is the flux density 
of a potential counterpart galaxy located a distance, r, away from 
the submillimetre position. 

A value of P < 0.05 is typically used as a threshold to se- 
lect robust radio IDs (e.g. Ivison et al. 2002, 2007) to submillimetre 
galaxies. We select all mid-infrared sources found within our search 
radii as potential counterparts and use the P statistic to distinguish 
between competing IDs when more than one galaxy is found. We 
note that almost all SMGs listed in Table[8]have mid-infrared coun- 
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1.0 




850 |a,m flux density (mJy) 850 flux density (mJy) 



Figure 3. Fractional completeness as a function of 850-/^m flux density for each of the five fields. The curves are fits to the data (see text for details). From 
left-hand panel shows the completeness for RX J163308 (squares), RX J121803 (diamonds) and RX J094144 (triangles) while the right-hand panel shows 
the completeness for RX J124913 (asterisks) and RX J005734 (diamonds). 




2 5 10 

850 ytim flux density (mJy) 



Figure 4. Integral number counts plotted in 1 mJy bins. The QSOs have 
been removed from the counts. The dotted line is a best-fit Schechter func- 
tion to the SHADES counts as described in Coppin et al. (2006) while the 
dashed line is the result of multiplying this fit by 4 and is shown for refer- 



terparts for which the computed P value is < 0.05 in at least one 
waveband, and can therefore be considered robust based on this 
criterion. 

The 4 QSOs detected at submillimetre wavelengths have coun- 
terparts in all 3 Spitzer bands. In addition, we find counterparts in at 
least one waveband for ^ 60 per cent of the SMGs. Postage stamp 
images of the mid-infrared counterparts are shown in Figs lAll and 
IA2l which can be found in Appendix A. 



7 MID-INFRARED COLOURS 

Mid-infrared colours provide a powerful means of distinguishing 
between galaxies which contain AGN and those that do not (e.g. 
Ivison et al. 2004, 2007; Egami et al. 2004). Luminous galaxies 
at 1 < z < 3 hosting AGN have power-law mid-infrared spectra 
whereas those without AGN are dominated by starlight in the IRAC 
bands and dust/PAH emission at 24 fim. 

In Fig.|5]we plot the ratio of 24 and 8 ^m flux density against 
the ratio of 8 and 4.5 flux density for those objects detected 
in at least two Spitzer bands. We superimpose redshift tracks for 
the local ULIRGs Mrk 231 (AGN-like SED) and Arp 220 and 
NGC 6240 (starburst-like SEDs). The QSOs and SMGs clearly in- 
habit a different regions of this diagram. As expected the QSOs 
have colours more similar to Mrk 231 than to the starburst galaxies. 
Three QSOs have remarkably similar mid-infrared colours given 
that they were selected on their X-ray properties. The other QSO, 
RX J 163308, has significantly different IRAC colours but we note 
from the postage stamp image (Fig. lAlb that it has a close 'com- 
panion' galaxy which is blended with the QSO in our extracted 
catalogue and will thus contaminate the colours. While the other 
three QSOs have 8.0/4.5 /im flux density ratios roughly consis- 
tent with the Mrk 231 template at their known redshifts we note 
that their 24/8.0 fim flux density ratios appear to be lower than that 
of Mrk 23 1 by a factor of ~ 2 which might be related to the rel- 
ative strength of emission/absorption features. All but one of the 
SMGs on the other hand have mid-infrared colours consistent with 
Arp 220 and/or NGC 6240 at 1 < 2 < 3. The one exception has 
colours consistent with Arp 220 at 0.1 < z < 0.5 but its faintness 
and non-detection at 24 jim are more consistent with a Mrk 231- 
like object at 1 < z < 3. It is thus the only SMG in our sample 
for which mid-infrared colours indicate significant emission from 
an AGN. We thus find no strong evidence for buried AGN in the 
SMGs. Having said this, we note that deep X-ray (e.g. Vignati et 
al. 1999; Netzer et al. 2005) and infrared spectroscopy (e.g. Lutz et 
al. 2003; Armus et al. 2006) observations of NGC 6240 show that 
this galaxy does in fact host a highly obscured AGN but that the 
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Table 8. Candidate Spitzer IRAC and MIPS counterparts to the submillimetre galaxies. Column 2 gives the ID number of the submillimetre galaxy as listed 
in Table|4]or Table|5] we searched for counterparts to all 450-/im detected objects that were not detected at 850 ^ra. Column 4 gives the catalogue number 
for each ID; for reference these numbers are also shown on Figs. lAll and lA2l Column 8 gives the angular' offset between the SCUBA position and the Spitzer 
position. Column 9 gives the connected P statistic (see text). 



Field 




ID 


A 

(^m) 


Cat* 


Sx 

(MJy) 


RA 
(J2000.0) 




Dec 
(J2000.0) 




A0 
(arcsec) 


P 


Note 


RXJ005734.78 


- 272827.4 


850#1 


4.5 


272 


476 ± 15 


QQh 57m 34s 


93 


—27° 


28' 


27' 


.6 


3.69 


0.005 


QSO 






850#1 


8.0 


203 


1534 ± 31 


QQh 57m 34s 


94 


—27° 


28' 


27' 


.6 


3.67 


0.003 


QSO 






850#1 


24.0 


10 


4107 ± 111 


00*" 57™ 34= 


90 


—27° 


28' 


27' 


.3 


4.12 


0.002 


QSO 






450#2 


4.5 


231 


26.5 ±3.6 


00^ 57™ 33= 


06 


—27° 


28' 


39' 


.7 


0.97 


0.009 








450#2 


8.0 


180 


18.4 ±4.4 


OQ^ 57™ 33= 


05 


—27° 


28' 


39' 


.6 


1.19 


0.011 




RXJ094144.51 


+ 385434.8 


850#1 


4.5 


242 


713 ± 19 


Qgh 4]^m 44s 


69 


+38° 


54' 


39' 


.8 


3.38 


0.003 


QSO 






850#1 


8.0 


180 


2147 ± 36 


Qgh 4]^m 44s 


70 


+38° 


54' 


39' 


.9 


3.41 


0.001 


QSO 






850#1 


24.0 


1 


5254 ± 182 


Ogh 4im 44s 


66 


+38° 


54' 


39' 


.5 


3.30 


0.001 


QSO 






850#2 


4.5 


191 


87.9 ±6.6 


og*" 41™ 46= 


34 


+38° 


53' 


56' 


.5 


2.80 


0.013 








850#2 


8.0 


155 


65.7 ±6.9 


ogh 41™ 46= 


33 


+38° 


53' 


55' 


.6 


2.45 


0.008 








850#2 


24.0 


6 


1041 ± 166 


ogh 41™ 46= 


38 


+38° 


53' 


57' 


.0 


2.89 


0.004 








850#3 


4.5 


339 


74.7 ±6.0 


ogh 41™ 42= 


66 


+38° 


55' 


36' 


.8 


2.53 


0.012 








850#3 


8.0 


241 


63.9 ±6.9 


ogh 41™ 42= 


68 


+38° 


55' 


36' 


.9 


2.50 


0.009 








450#3 


4.5 


205 


19.0 ± 3.1 


ogh 41™ 47= 


00 


+38° 


54' 


11' 


.2 


4.90 


0.124 








450#3 


4.5 


206 


13.1 ± 2.6 


ogh 41™ 46= 


39 


+38° 


54' 


13' 


.9 


4.11 


0.131 








450#3 


8.0 


159 


24.9 ±4.8 


ogh 41™ 46= 


95 


+38° 


54' 


10' 


.9 


4.69 


0.055 








450#3 


8.0 


160 


22.3 ±4.7 


ogh 41™ 46= 


39 


+38° 


54' 


13' 


.9 


4.11 


0.052 




RXJ121803.82 


+ 470854.6 


850#2 


4.5 


223 


47.8 ±4.9 


12h 18™ 07= 


57 


+47° 


08' 


45' 


.2 


0.92 


0.008 








850#2 


8.0 


175 


35.9 ±5.3 


12h 18™ 07= 


57 


+47° 


08' 


45' 


.0 


0.90 


0.004 








850#2 


24.0 


6 


476 ± 104 


12h 18™ 07= 


60 


+47° 


08' 


46' 


.8 


1.76 


0.001 








850#4 


4.5 


256 


4.5 ± 1.5 


12h 18™ 02= 


54 


+47° 


08' 


28' 


.3 


4.07 


0.317 








850#4 


8.0 


188 


16.7 ±4.0 


12h 18™ 02= 


67 


+47° 


08' 


27' 


.2 


2.48 


0.038 








850#5 


24.0 


5 


3200 ± 160 


12h 18™ 04= 


46 


+47° 


08' 


50' 


.6 


8.87 


0.006 


likely false ID 


RXJ124913.86 


- 055906.2 


850#1 


4.5 


261 


3671 ± 42 


12h 49™ 13= 


88 


-05° 


59' 


19' 


.0 


1.71 


< 0.001 


QSO 






850#1 


8.0 


200 


12764 ± 88 


12h 49™ 13= 


90 


-05° 


59' 


18' 


.9 


1.46 


< 0.001 


QSO 






850#1 


24.0 


5 


26311 ± 125 


12h 49™ 13= 


86 


-05° 


59' 


18' 


.4 


2.13 


0.001 


QSO 






850#4 


4.5 


342 


18.1 ± 3.0 


12h 49™ 12= 


56 


-05° 


58' 


34' 


.9 


3.18 


0.093 








850#4 


4.5 


343 


79.1 ± 6.2 


12h 49™ 12= 


26 


-05° 


58' 


30' 


.6 


3.34 


0.018 








850#4 


8.0 


265 


55.9 ± 7.5 


12h 49™ 12= 


26 


-05° 


58' 


31' 


.2 


2.82 


0.015 








850#5 


4.5 


324 


25.5 ± 3.6 


12h 49™ 11= 


39 


-05° 


59' 


57' 


.9 


2.89 


0.061 








850#5 


8.0 


236 


33.7 ± 6.6 


12h 49™ 11= 


44 


-05° 


59' 


58' 


.4 


2.00 


0.023 








450#2 


4.5 


404 


2.2 ± 1.1 


12h 49™ 10= 


95 


-05° 


58' 


34' 


.7 


3.06 


0.194 




RXJ163308.57 


+ 570258.7 


850#1 


4.5 


142 


21.3 ± 3.3 


16h 32™ 59= 


74 


+57° 


03' 


21' 


.5 


3.48 


0.076 








850#1 


8.0 


134 


30.0 ± 5.3 


16h 32™ 59= 


74 


+57° 


03' 


21' 


.3 


3.29 


0.037 








850#1 


24.0 


17 


275 ± 55 


16h 32™ 59= 


65 


+57° 


03' 


20' 


.8 


3.47 


0.017 








850#2 


4.5 


242 


83.3 ± 6.4 


16h 33™ 08= 


55 


+57° 


02' 


54' 


.6 


4.67 


0.038 


QSO 






850#2 


8.0 


226 


153 ± 10 


16h 33™ 08= 


52 


+57° 


02' 


54' 


.5 


4.86 


0.023 


QSO 






850#2 


24.0 


18 


705 ± 40 


16h 33™ 08= 


54 


+57° 


02' 


53' 


.0 


6.18 


0.004 


QSO 






850#5 


4.5 


174 


19.1 ±3.1 


16h 33™ 03= 


17 


+57° 


02' 


58. 


"6 


0.13 


< 0.001 








850#5 


8.0 


168 


17.0 ±4.5 


16h 33™ 03= 


18 


+57° 


02' 


58' 


.9 


0.16 


< 0.001 








850#5 


24.0 


19 


256 ± 32 


16h 33™ 03= 


09 


+57° 


02' 


57' 


.6 


1.27 


0.002 








850#6 


4.5 


290 


2.9 ± 1.3 


16h 33™ 12= 


45 


+57° 


02' 


34' 


.9 


3.20 


0.172 








850#7 
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Figure 5. Colour-colour plot for those sources detected in two or three of the observed Spitzer bands (24, 8.0 and 4.5 ^m). The QSOs are plotted as diamonds 
while the submiUimetre-selected galaxies discovered in their fields are shown as squares. Lines show the redshifted colour tracks of three local ultraluminous 
infrared galaxies with (Mrk 231) and without (Arp 220 and NGC 6240) obvious signatures of powerful AGNs at mid-infrared wavelengths. Redshifts of 1.0, 
2.0 and 3.0 are marked on each of these tracks. 



energetics are dominated by the starburst. This leaves open the pos- 
sibiUty that many of the SMGs contain buried AGN which are not 
powerful enough to dominate over the starburst in the mid-infrared. 



8 DISCUSSION AND CONCLUSIONS 

We have shown that the number counts of SMGs with flux densi- 
ties ~ 2 — 4 njy in the ~ 2 arcmin diameter fields surrounding our 
X-ray absorbed QSOs is higher than that found for a blank field. 
Most of these galaxies have mid-infrared colours very similar to 
the local ULIRG Arp 220 if it were placed at 1 < 2 < 3. Fitting 
the millimetre to far-infrared spectrum of Arp 220 with a greybody 
and then integrating to determine the far-infrared luminosity we 
find LpiR = (1.3 ± 0.1) X lO^'^ which corresponds to a star- 
formation rate (SFR) of 224± 17 M0 yr^^ after applying the usual 
scaling-law (Kemricutt 1998). Using the same greybody fit and ap- 
plying the appropriate k-correction leads to predicted 850 jjm flux 
densities of 1.73, 1.44 and 1.42 mjy for Arp 220 at z=l, 2 and 3 
respectively. The SMGs in our sample have 850 /im flux densities 
ranging from about 3—10 mJy implying that their SFRs range be- 
tween 400 — 1300 Mq yr^^ if Arp 220 is indeed a good analogue. 
Therefore, if the starburst can be sustained for a few hundred mil- 
lion years or so then a stellar mass equivalent to that of the bulge of 
a large galaxy will be assembled. 

It is natural to speculate that an over-dense region of space at 
high redshift containing galaxies that are capable of building large, 
homogeneous populations of stars wiU evolve over cosmic time to 



become clusters of galaxies containing a population of massive el- 
Uptical galaxies. Since it is well established that elliptical galax- 
ies contain dormant SMBHs (e.g. Magorrian et al. 1998), and that 
these must have assembled at high redshifts where the accretion 
luminosity density peaks, we expect that our SMGs should con- 
tain growing black holes, i.e. AGN. How the major growth phase 
of the black hole coincides temporally with the luminous starburst 
that builds up the stellar mass is a current topic of debate. Indeed, 
models predict a range of possibilities (e.g. Archibald et al 2002; 
Di Matteo, Springel & Hernquist 2005) but there must be a link 
between the formation of the stars in the bulge and the SMBH be- 
cause of the tight correlation between their masses found for local 
galaxies (Ferrarese & Merritt 2000). After a number of early nuU 
results using moderate depth X-ray data (e.g. Fabian et al. 2000; 
Waskett et al. 2003) deep X-ray (Alexander ct al. 2005a,b, 2008) 
and mid-infrared spectroscopy (Menendez-Delmestre et al. 2007; 
Pope et al. 2008) observations show that SMGs found in the field 
often do contain AGN but that the masses of their black holes may 
be much smaller than those found in high redshift QSOs, and con- 
sequently that the bulk of the far-infrared emission is due to the 
starburst. Moreover, estimates of the host galaxy masses suggest 
that the growth of the black hole lags the growth of the spheroid 
(Borys et al. 2005; Alexander et al. 2008; see also Coppin et al. 
2009). While this may be the case for SMGs, analysis of massive 
galaxies selected on other criteria find evidence to the contrary, sug- 
gesting that more than one scenario is possible for the formation of 
a massive galaxy and its SMBH. For example, studies of powerful 
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radio sources (McLure et al. 2006) and luminous QSOs (Peng et 
al. 2006) suggest that the black-hole mass to bulge luminosity ratio 
increases towards high redshifts thus indicating that the black hole 
reaches its critical mass before the galaxy bulge is fully formed. 

The mid-infrared colours of our SMGs also argue against the 
presence of a powerful AGN in these objects unless they are very 
highly obscured. They are thus either intrinsically different objects 
to the QSO or they are caught at a different time in their evolution. 
Page et al. (2010, in prep.) argue that the X-ray absorption in the 
QSOs has its origin in an ionized absorber associated with a power- 
ful wind capable of terminating the star formation. Such feedback 
is a popular ingredient of models since it leads naturally to the rela- 
tionship between SMBH mass and spheroid mass discussed above 
(e.g. Fabian 1999). If the SMGs in the same field evolve to form 
a cluster of ellipticals then they would presumably go through a 
similar evolutionary sequence to the QSO but are caught at an ear- 
lier stage of their life cycle at 1 < 2; < 3. Estimates of the stellar 
masses would be useful to put further constraints on their evolu- 
tionary status. To this end we are currently fitting models to optical 
through infrared data to determine photometric redshifts and rest- 
frame K-h&nd luminosities; an analysis of the RXJ 094144 field is 
currently in progress (Carrera et al. 2010, in prep.). 

One difference worthy of note between the findings of this 
paper and those for high-redshift (typically 2; > 3) radio galax- 
ies (HzRGs) concerns very bright SMGs in their fields. Two of 
the seven HzRG fields imaged with SCUBA contain SMGs with 
Sgso > 20 mjy (Stevens et al. 2003) whereas the brightest ob- 
ject found around the QSOs has S'sso ~ 10 mJy. Although the 
samples are too small to give a statistically significant result it ap- 
pears that the SMGs detected around the lower redshift, less pow- 
erful, radio-quiet AGN have lower Sgso- A direct comparison is 
not informative because the sources were extracted with different 
signal-to-noise criteria and by different methods but comparing flux 
densities significantly above the catalogue cut-offs yields 8 SMGs 
with Sgso > 6 mJy in the 7 HzRG fields while only 2 are detected 
around the 5 QSOs. 

Given the effectively flat relationship between Sgso and red- 
shift for 1 < jz < 10 this difference in flux density leads to a 
proportional difference in star-formation rates (all other things be- 
ing equal). There is thus reasonably good evidence for a difference 
in star-formation activity in individual galaxies in these fields in- 
dicating that the galaxy formation process might be dependent on 
the magnitude of the density peak and/or its redshift. For exam- 
ple, in the case of HzRGs it was found that the detection rate in- 
creased dramatically at 2 > 2.5 (Archibald et al. 2001; Reuland 
et al. 2004) which is significantly higher than the average redshift 
of our QSO sample. Another possibility is that the radio proper- 
ties of the central AGN are in some way linked to environmental 
density (e.g. Kauffmann, Heckman & Best 2008). Indeed, recent 
work by Falder et al. (2010) finds evidence for larger over-densities 
of galaxies around radio-loud objects (based on work carried out 
at 3.6 /im). Perhaps radio synchrotron luminosity is enhanced in 
regions where the inter-galactic medium is more dense. 

An obvious next step will be to conduct a statistically mean- 
ingful survey of fields around AGN over the full range of red- 
shifts in order to track star-formation activity as a function of cos- 
mic time. The next generation bolometer camera on the JCMT 
(SCUBA-2) and SPIRE/PACS on-board the Herschel Space Obser- 
vatory are instruments with good enough sensitivity to make such a 
campaign practical. It would also allow a bigger area to be mapped 
around each AGN, better matching the size of proto-cluster regions 
predicted by numerical simulations. 
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APPENDIX A: NOTES AND POSTAGE STAMP IMAGES 

In this section we present and discuss postage stamp images of 
the QSOs and SMGs. For the RX J094144 field we compare the 
mid-infrared counterparts to the SMGs with those identified pre- 
viously at optical and near-infrared wavelengths. We also discuss 
cases where likely mid-infrared counterparts were not catalogued 
because they fell just outside our formal search radii, and one case 
where an obvious SMG was not found by our extraction algorithm; 
these sources were not included in the analysis described in the 
main part of this paper. 

Al The RX J094144.51+385434.8 field 

The SCUBA data for this field were published previously by 
Stevens et al. (2004). In the present paper we have treated the 
data slightly differently in the final reduction stages (rebinning and 
smoothing). We also apply a formal source extraction algorithm in 
this work resulting in only 2 sources at 850 /im and 4 at 450 nm 
whereas Stevens et al. (2004) simply identified peaks on the maps 
above a 3-a signal-to-noise threshold, identifying 6 sources. As al- 
ready noted in Section 143] the source extraction algorithm fails to 
detect sources that are blended which has undoubtedly affected the 
850 fim extraction. At 450 nm we can see that sources 450#1 and 



450#3 appear more elongated than the other two objects hinting 
that they may be resolved. It was these objects that Stevens et al. 
(2004) identified as blends. 

Inspection of the postage stamps shows that source 
850#2/450#2 is coincident with a bright mid-infrared source de- 
tected in all 3 Spitzer channels. This object was one of a pair of 
sources detected in the R- and Ji'-bands by Stevens et al. (2004). 
The new Spitzer data show that the correct identification is the ob- 
ject with the redder R — K colour since this object becomes pro- 
gressively more prominent with increasing wavelength compared 
to its neighbour and is detected at 24 pm. 

Source 450#3 falls in a crowded field. We find 2 possible 
identifications in the IRAC channels. Of these source 206/160 lies 
closer to the 450 pm source centroid and is also coincident with 
faint (formally undetected) emission at 24 pm. Source 206/160 
was not one of the possible identifications discussed by Stevens 
et al. (2004) since it is very faint/undetected in the J-f-band. We 
conclude that the likely identification for this SMG is the very red 
object 206/160. 

Source 450#4 was matched with an extremely red object 
(R ~ K = 5.8) by Stevens et al. (2004). This source is not for- 
mally recovered in this work because it lies outside our 5 arcsec 
IRAC search radius but is visible in the postage stamp to the west 
of the 450 ^m centroid as a bright source at 4.5 and 8.0 /im. Given 
its red colours we consider this object as a likely identification. 
Measured flux densities are 57.0 ± 5.3 and 48.5 ±6.1 /iJy at 4.5 
and 8.0 respectively giving a colour similar to Arp 220 at the 
redshift of the QSO. The 4.5-/im centroid is 6.6 arcsec from the ex- 
tracted 450 pm position at 09*' 41™ 43". 61, -^38° 54' 24". 5 while 
the 8.0-/im centroid is 6.4 arcsec from the extracted 450 pm posi- 
tion at 09'' 41™ 43=. 64, +38° 54' 24". 6. The calculated P-statistic 
is 0.089 and 0.064 at 4.5 and 8.0 pm respectively. 



A2 The RX J121803.82+470854.6 field 

In Section |431 we noted that the QSO in this field was not found by 
our extraction algorithm but was detected by Page et al. (2001). The 
latter observations used SCUBA in 'photometry' mode where the 
central bolometer of each array is centred on the target but no image 
is made. Inspection of Fig.[T]shows that while there is submillime- 
tre emission at the location of the QSO (located below source 5 
on the 850 pm image) the surrounding structure is complex. It is 
thus possible that the source extraction algorithm failed to locate a 
significant peak because of blending. It is also likely that the pho- 
tometry observation of this object reported by Page et al. (2001) is 
contaminated by emission from the SMGs close to the QSO. 

In Table[8]we noted that object 850#5 has a likely false iden- 
tification at 24 pm. This can also be attributed to the crowded field 
close to the QSO since it is the QSO itself that has fallen within the 
search radius at 24 pm while it was not picked up in the smaller 
search radii adopted for the IRAC channels. 



A3 The RX J124913.86 055906.2 field 

Source 850#2 has no listed identification in Table[8]but the postage 
stamp images show bright 4.5 and 8.0 ^m sources close to the 
850 pm centroid. The closest 4.5 pm source is at 12*^ 49™ 15M9, 
-05° 58' 40". 9 which at 5.3 arcsec from the SMG is only just 
outside our formal search radius. This object has a flux density of 
17.7±3.0 piy and is coincident with an 8.0 ^m source although the 
latter is not detected by SEXTRACTOR. Inspection of the 8.0 pm 
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image reveals that the diffraction spike from the bright QSO at the 
field centre falls on top of this object probably explaining the non- 
detection. Although the P-statistic for the 4.5 /.im source is high at 
0.259 we consider this object as a likely identification because of 
its detection in both IRAC bands. 

Sources 850#4 and 850#4 have an obvious identification in 
the IRAC bands as well as faint (formally undetected) emission at 
24 jim. 

A4 The RX J163308.57+570258.7 field 

Source 850#4 has no listed identification in Table[8]but the postage 
stamp images show a bright 4.5 and 8.0 /im source close to the 
850 /im centroid along with faint (formally undetected) emission 
at 24 /im. The 4.5 /xm source is located at 16'' 33"" 10M5, 
+57° 02' 10". 2 only 5.2 arcsec from the SMG, has a flux den- 
sity of 12.9 ± 2.6 and a P-statistic of 0.266. The 8.0 /xm source is 
located at 16*' 33™ 10". 08, +57° 02' 09". 3 only 5.9 arcsec from 
the SMG, has a flux density of 17.6 ±4.5 and a P-statistic of 0. 198. 
This object has IRAC colours similar to Arp 220 at the redshift of 
the QSO. The proximity of this object to the SMG together with its 
detection in in both IRAC bands make it a likely identification. 

Inspection of the 850 /xm map for this field (Fig. [T) 
shows submillimetre emission north of the QSO at approximately 
le*" 33™ 05=. 75, +57° 04' 11". 7. There is a negative region to 
the east of this position which is consistent with a negative image 
of the source located at one of the off-positions. The western off- 
position would fall off the map. We believe that this is a real SMG 
(probably two sources) which was missed by the extraction algo- 
rithm because of its proximity to the edge of the map. Inspection of 
the postage stamp images reinforces this belief since we find two 
24 fj,m sources coincident with the SMG position and at the same 
approximate position angle as the 850 /im emission. The IRAC data 
reveal 3 and 5 sources at 4.5 and 8.0 jim respectively. The southern 
24/850 /im source has two possible IRAC identifications. Of these, 
source 260/232 is the brightest and reddest object; its 8.0/4.5 /xm 
flux density ratio is 1.96 ± 0.22. If source 10 is the 24 /im counter- 
part then the ratio of flux densities at 24 and 8.0 fjm is 4.60 ± 0.99 
which would place this object on the Mrk 231 curve shown in Fig|5] 
indicating that it hosts a buried AGN, possibly at a lower redshift 
than RX J163308. The northern 24/850 /xm source has 1 counter- 
part at 4.5 ijm and 3 at 8.0 /xm although the 4.5 fjm emission is 
extended in the same manner as that at 8.0 ^m implying that this 
SMG could have multiple mid-infrared counterparts. Inspection of 
the flux densities presented in Table[8]shows that the IRAC colours 
are similar to those of Arp 220 at the redshift of the QSO so this 
source is unlikely to host an AGN. Since this system is one of the 
most interesting in our sample we include it in Table [8] as source 
850#7 although it is not included in the 850 ^m number counts 
analysis because it was not found by the extraction algorithm. 

This paper has been typeset from a TgX/ KTeX file prepared by the 
author. 
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Figure Al. Postage stamps of 25 X 25 arcsec cut out at the extracted right ascension and decUnation of the 850 fira sources from Table|4] The panels show 
(left) MIPS 24 /im contours (black) overlaid on the IRAC 4.5 /^m data (grey-scale and white contours) and (right) SCUBA 850 /tm contours (black) overlaid 
on the IRAC 8.0 /tm data (grey-scale and white contours). Contour levels at 850 ^m start at ~ 2 cr and increase in 1 cr steps while the Spitzer contours are 
chosen to highlight faint sources that might not have been formally extracted. The numbering corresponds to the catalogue numbers given in Table[8] 
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Figure Al. - continued 
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Figure A2. Postage stamps of 25 X 25 arcsec cut out at the extracted right ascension and decUnation of the 450 fira sources from Table |5] The panels show 
(left) MIPS 24 /im contours (black) overlaid on the IRAC 4.5 /^m data (grey-scale and white contours) and (right) SCUBA 450 /tm contours (black) overlaid 
on the IRAC 8.0 /jm data (grey-scale and white contours). Contour levels at 450 ^tm start at ~ 2 cr and increase in ~ 1 cr steps while the Spitzer contours are 
chosen to highlight faint sources that might not have been formally extracted. The numbering corresponds to the catalogue numbers given in Table[8] 



